The application of slim disk models to ULX: 
the case of M33 X-8 
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^ ■ Abstract 

A comparative XMM-Newton archival data spectral study of the ultraluminous X- 



ray source (ULX) M33 X-8 has been performed by using both the standard disk 



model and the newly developed slim disk models. The results of this analysis tend 
to confirm the hypothesis that M33 X-8 is an X-ray binary with a stellar-mass black 
hole accreting at super-Eddington rate. 
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Since the early observations with the Einstein satellite, many ultraluminous 
X-ray sources (ULX) have been discovered (see, for example, Fabbiano 1989, 
Colbert & Mushotzky 1999, Makishima et al. 2000, Colbert & Ptak 2002, 
Foschini et al. 2002a, Miller & Colbert 2004). ULX are point-like sources 
located sufficiently far from the centre of the host galaxy and with X-ray 
luminosity of about 10 39-40 erg/s in the 0.5 — 10 keV energy band. According 
to their X-ray spectrum, ULX can be broadly divided into three classes: 
Type I: a thermal component, usually modelled with a multicolour disk (MCD) 
with an inner disk temperature T in « 1.0—1.5 keV, plus a power-law-like excess 
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at high energies. 

Type II: similar to Type I, but with an inner disk temperature T in 0.2 keV; 
these ULX are candidates to host intermediate-mass black holes (e.g. Miller 
et al. 2003, 2004); 

Type III: a simple power-law model, with a photon index T ~ 2; these ULX 
could be either background objects (e.g. NGC4698 X-l, Foschini et al. 2002b; 
NGC4168 X-l, Masetti et al. 2003; see Foschini et al. 2002a for a discussion on 
the probability to find background AGN in ULX surveys) or stellar-mass black 
holes in an anomalous very high state (Comptonization-dominated, Kubota 
et al. 2002). 

It is worth mentioning that some ULX show state transitions, implying also a 
transition in the types of the present classification (e.g. Kubota et al. 2001a). 

The model mainly used to fit the thermal component is based on the standard 
optically thick and geometrically thin accretion disk by Shakura & Sunyaev 
(1973). The multicolour disk version developed by Mitsuda et al. (1984) is 
widely applied in the fit of ULX spectra (diskbb model in xspec). However, 
this model, together with the power-law model used to fit the hard X-ray 
excess, presents some problems addressed by several authors (Merloni et al. 
2000, Ebisawa et al. 2003, Wang et al. 2004). Wang et al. (2004) proposed a 
Comptonized version of the multicolour disk (CMCD) to take into account the 
feedback effects of the disk-corona system and applied this model to some well- 
known ULX. They found an apparent agreement between the results obtained 
with the MCD+PL and the CMCD model, despite the known limitations. 

On the other hand, the emission from a number of Galactic black hole candi- 
dates (e.g. GRO J1655-40, Kubota et al. 2001b, XTE J1550-564, Kubota et 
al. 2002) cannot be explained in terms of the canonical Ldisk oc relationship 
between the accretion disk bolometric luminosity and its inner temperature. 
This behaviour has been interpreted in terms of an anomalous state which 
can be represented with an optically and geometrically thick accretion disk, 
or "slim disk" (e.g. Watarai et al. 2000). Since this type of disk can support 
super-Eddington luminosities, the slim disk model was also proposed to fit 
some ULX spectra (e.g. Watarai et al. 2001, Kubota et al. 2002, Ebisawa et al. 
2003). If this model is correct, some ULX could be understood as stellar-mass 
black holes with supercritical accretion rates, thus avoiding the introduction 
of intermediate-mass black holes. 

In the present work, we focus on a single ULX - M33 X-8 (Type I) - and 
perform a comparative X-ray energy spectral study by using both the standard 
disk and the newly developed slim disk models. Specifically, for the latter, 
there are several codes available (e.g. Mineshige et al., 2000 without relativistic 
effects and electron scatterings; Watarai et al. 2000, without relativistic effects 
and a first order approximation for the electron scattering), but we adopted the 
version developed by Kawaguchi (2003), where the effects of electron scattering 
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and relativistic correction are included. Public archival data, obtained with 
XMM-Newton observations, are fit by these models in the present study. 

M33 X-8 can be considered an object between a well-known binary system 
(such as Cygnus X-l) and a typical ULX (with L x ~ 10 395 erg/s). Therefore, 
it could be a good laboratory to start this type of comparative study. Since 
the discovery of M33 X-8 with the Einstein satellite (Long et al. 1981), it was 
clear that the source presented peculiar properties (Lo.2-4 kcv ~ 10 39 erg/s, 
soft spectrum) and already Trinchieri et al. (1988) suggested that it could 
be a new type of X-ray binary system. Takano et al. (1994) on the basis of 
ASCA data, suggested that it could be a stellar mass BH (ps 1OM ), although 
they were skeptical on the possibility to find such a peculiar binary system so 
close to the M33 centre. Makishima et al. (2000) suggested that X-8 can be 
considered as a member of ULX population. 

A detailed analysis of a set of XMM-Newton observations of X-8 has been 
presented by Foschini et al. (2004). The statistically best-fit model is the 
MCD+PL, in agreement with the findings of other authors (e.g. Makishima 
et al. 2000, Takano et al. 1994), but several correction factors were needed to 
infer useful physical parameters. Specifically, for the case of a face-on disk, the 
inferred mass is (6.2 ± 0.4)M Q , which can increase up to (12 ± 1)M if the 
correction factors for Doppler boosting, gravitational focusing, and a 9 = 60° 
disk inclination, are taken into account. 

For the present work, we retrieved from the XMM-Newton public archive a set 
of three observations with M33 X-8 on axis (ObsID 0102640101, 0141980501, 
0141980801), performed on 4 August 2000, 22 January and 12 February 2003, 
respectively. The data were reduced and analyzed in a manner similar to that 
described in Foschini et al. (2004), except for the use of a more recent version of 
the software (XMM-SAS v . 6.1.0 and HEASof t v 6.0). Moreover, we analyze 
here for the first time the 2003 observations, now become public. The spectra 
were fit in the energy range 0.4 — 10 keV for PN and 0.5 — 10 keV for MOS1 
and MOS2 (in this case, the flux was then extrapolated down to 0.4 keV). 
The width of the energy range was selected to take into account the present 
status of the instrument calibration (Kirsch 2004). The results of the fits are 
presented in Table 1 and Figure 1. 

From a first look at the spectral fits, the model of Shimura & Takahara (1995) 
results in a varying black hole mass. The GRAD model is not able to fully 
account for the high-energy emission of M33 X-8, thus resulting in a clear 
excess for E > 5 keV. The same occurs with the MCD model alone (not 
shown), but by adding a power-law component it is possible to model these 
excesses (see the MCD+PL fit). The slim disk model SDK provides a fit 
comparable to the MCD+PL, with constant mass within the measurement 
errors. However, it should be noted that in order to fit all the data with the 
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Fig. 1. XMM-Newton-EPIC spectra of M33 X-8 (ObsID 0102640101) and fit residu- 
als to: (top-left) general relativistic accretion disk (grad model in xspec); (top-right) 
multicolour accretion disk (diskbb model) plus a power-law model; (bottom-left) 
slim disk model by Kawaguchi; (bottom-right) Shimura & Takahara's model. 

MCD+PL models for a constant BH mass, fine-tuning of some parameters is 
needed (spectral boosting factor and inner radius of the disk etc.). In addition, 
we also note that the luminosity of the PL component at 0.4 — 10 keV range 
dominates over the disk luminosity. This is not usual for Galactic BH, although 
two other ULX show similar behaviour (NGC55 X-7, Stobbart et al. 2004; 
NGC5204 X-l, Roberts et al. 2005). 



The results presented here are only preliminary results of a comparative study 
on the X-ray spectra of ULX: it is necessary still to improve the models and to 
extend a more detailed analysis on a larger sample of ULX. Particularly, the 
spectral behaviour vs time appears to be the key to understand the nature of 
ULX. 
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Table 1 

Results from the fit of the X-ray data. Models: General relativistic accretion disk 
(GRAD; grad model in xspec), multicolour accretion disk (MCD; diskbb model in 
xspec), power-law model (PL), Kawaguchi's (2003) slim disk model (SDK), Shimura 
& Takahara's (1995) model (ST) with a fixed to 0.1. Meaning of symbols: Ah 
absorption columns [10 21 cm -2 ], with the Galactic column Ah = 5.6 x 10 20 cm -2 ; 
mass M [Mq] (for the MCD+PL model, the mass is calculated from the diskbb 
model normalization); accretion rate M [LEdd/c 2 ]; viscosity parameter a; observed 
flux F in the 0.4 — 10 keV energy band [10 -11 erg cm" 2 s -1 ]; intrinsic luminosity 
L in the same energy band [10 39 erg/s], calculated for d = 795 kpc; photon index 
T; inner disk temperature T in [keV]. We considered only the case of face-on disk. A 
normalization constant between MOS1/MOS2 and PN in the range 0.83 — 0.91 has 
been used. The uncertainties in the parameters are at the 90% confidence level. 



Model 


Parameter 


0102640101 


0141980501 


0141980801 


GRAD 


A H 


0.56 (fixed) 


0.56 (fixed) 


0.56 (fixed) 




M 


5.7 ±0.2 


5.4 ±0.3 


6 5 +0 - 3 

D.O_q 2 




M 


18.9 ±0.3 


21.1 ±0.5 


16.5 ±0.3 




F/L 


1.8/1.5 


2.0/1.6 


1.5/1.2 




X 2 /dof 


1.09/763 


1.12/361 


1.55/650 


SDK 


N H 


1.26 ±0.06 


1 q+0.2 
1 -°-0.1 


1.00 ±0.07 




M 


6 4 +0 - 6 


7+11 
'-2 


8±1 




M 


31 ±3 


QO + 54 


2211 




a 


< 0.005 


< 0.01 


< 0.01 




F/L 


1.8/1.7 


2.1/1.9 


1.7/1.5 




X 2 /dof 


1.00/761 


0.99/359 


1.15/648 


ST 


N H 


0.78 ± 0.06 


0.9 ±0.1 


0.56 (fixed) 




M 


10 l +ai 

lu - i -0.3 


10-118:? 


8.0 ±0.3 




M 


11 2 +a3 


12.818:? 


11.8 ±0.4 




F/L 


1.8/1.6 


2.1/1.8 


1.7/1.4 




X 2 /dof 


1.02/762 


1.03/360 


1.34/650 


MCD+PL 


N u 


2.2 ±0.4 


2 7+ 1 - 1 


9 -1+0.3 




M 


c c+0.9 
°-°-0.6 


O O + 0.9 

d - 8 -0.8 


o q+42.7 
6.S_ 27 






i 91 +0.08 
1 - zi -0.09 


1 5 +0 - 2 
± -°-0.3 


9 +0 - 7 
u - y -0.3 




r 


2 6 +0 - 3 

Z.D_q 2 


9 7 +0.9 
z -'-0.5 


2.2 ±0.2 




F/L 


1.8/2.1 


2.1/2.5 


1.7/1.9 




X 2 /dof 


1.00/760 


1.00/358 


1.09/647 
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